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Abstract 

We consider a single Minkowski brane sandwiched in between two copies of 
anti-de Sitter space. We allow the bulk Planck mass and cosmological constant 
to differ on either side of the brane. Linearised perturbations about this back- 
ground reveal that gravity can be modified in the infra-red. At intermediate 
scales, the braneworld propagator mimics four-dimensional GR in that it has 
the correct momentum dependance. However it has the wrong tensor struc- 
ture. Beyond a source dependant scale, we show that quadratic brane bending 
contributions become important, and conspire to correct the tensor structure 
of the propagator. We argue that even higher order terms can consistently be 
ignored up to very high energies, and suggest that there is no problem with 
strong coupling. We also consider scalar and vector perturbations in the bulk, 
checking for scalar ghosts. 
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1 Introduction 



How much do we believe in Einstein's General Theory of Relativity? We certainly 
don't believe it at energies beyond the Planck scale, where a quantum theory is 
required. In fact, experimental tests of gravity have only ever been conducted at 
much lower scales. For example, Newton's law has been tested down to about 0.2 
mm, or equivalently up to energies of around 10~^ eV. Solar system tests of GR are 
only valid outside of the Schwarzschild radius of the Sun. 

Ultra-violet modifications of gravity beyond 10~^ eV occur in all higher dimen- 
sional theories, such as String Theory. Recently, the idea that gravity is also modified 
in the infra-red has received much interest [U I2I El IH El El El • This is entirely plau- 
sible. Again, Newton's law has only been tested up to distances of around 10^^ cm, 
and solar system tests of GR do not extend much beyond the size of the solar system 
itself. 

Infra-red modifications of gravity have some attractive features. They can be 
used to tackle the cosmological constant problem jHl E] by weakening gravity at low 
energies. To understand this in a simple way, suppose we have Newton's constant G 
at intermediate scales, and Gir <^ G in the infra-red. Even for a large cosmological 
constant A, the contribution to the Hubble rate is considerably suppressed at large 
values of the scale factor 

H"^ ~ GirA < GA (1) 

However, that is not to say that these models prohibit late time inflation. On the 
contrary, these models generically give rise to Friedmann equations that are non- 
linear in H^. These permit self accelerating solutions that can be tuned to the current 
Hubble rate [IHl HH IH ESI El CSl • 

Unfortnately, these models are not without their problems. The most common 
of these is the vDVZ discontinuity that appears in the graviton propagator at all 
scales |lti| I17j. This represents a deviation from GR that will give the wrong predic- 
tion for light bending around the Sun. In some cases, it has been argued that this 
discontinuity reflects a breakdown in the linearised analysis jlHl EH] , as opposed to 
a genuine sickness in the theory. Some theories predict ghost states, and therefore 
violate unitarity [20] • A sickness as bad as this usually leads to the death of a theory. 
We should emphasise, however, that some models are indeed ghost-free [H Ej • 

In we considered asymmetric branes. By asymmetric, we mean that the 
gravitational parameters of the theory can differ on either side of the brane. There 
are at least two ways in which this asymmetry might arise. Firstly, suppose we 
have some sort of wine bottle shaped compactification down to 5 dimensions. In 
the effective theory, the Planck scale at the fat end of the bottle will be less than 
that at the thin end. Secondly, in |22l EH] we showed how to construct a domain 
wall living entirely on the brane. In some cases the Planck scale on the brane 
differed on either side of the domain wall. Asymmetric branes admit self accelerating 
solutions 121]. It is natural to ask, therefore, if asymmetric brane models exhibit 
infra-red modifications of gravity. 

In this paper, we will consider a positive tension brane. To the left of the brane, the 
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bulk is described by the Einstein-Hilbert action with Planck mass, M^^ and negative 
cosmological constant, = ~6//|. To the right, we have Planck mass, Mr, and 
cosmological constant, A/j = —Q/lj^. In general, we will take 3> Mji and 1^ -C Ir. 

We will be particularly interested in Minkowski branes. Because of the asymmetry, 
we are allowed three classes of solution. We have (i) a Randall-Sundrum like solution 
where the warp factor in the metric decays away from the brane on both sides |24], 
(ii) an "inverse" RS solution where the warp factor grows on both sides, and (iii) a 
mixed solution where the warp factor decays on one side and grows on the other. 

If we place a small matter source on the brane, we can investigate linearised 
perturbations about the background solutions. For the RS solution and the mixed 
solution we see the introduction of a new length scale 



At intermediate energies (1/r <C p -C the braneworld propagator is propor- 

tional to and will reproduce Newton's law. At lower energies {1/Ir -C p -C 1/r), 
the momentum dependance changes, so that Newton's law is modified in the infra-red. 
However, as in the DGP model P], the propagator suffers from a vDVZ discontinuity 
at intermediate scales, and will predict the wrong results for light bending around 
the Sun. 

It turns out that the vDVZ discontinuity can be removed by including brane 
bending terms to quadratic order. For a source of mass m, we discover yet another 
scale 



For p^: p <^ the quadratic terms dominate over the linear ones. The re- 

sult is that the braneworld propagator agrees exactly with four- dimensional General 
Relativity. We might be worried about the break down of the linearised analysis at 
this scale. In the DGP model, this breakdown has been linked to a strong coupling 
problem [23 12H1 123 12H] • However, we will argue that there is nothing to worry about. 
This is because even though quadratic terms become important when p ~ p*, even 
higher order terms are only important at much higher energies. It seems that we can 
indeed consistently modify four-dimensional GR in the infra-red. 

The rest of this paper is organised as follows. In section |2l we describe our set-up 
in detail, before deriving background solutions in section El In section 0] we derive 
the linearised equations of motion and solve them. We analyse the solutions and 
establish the existance of the long distance scale r, aswell as the vDVZ discontinuity. 
In section |SJ we consider scalar and vector perturbations in the bulk. Vector pertur- 
bations only exist for the inverse RS solution. Scalar perturbations exist whenever 
we include the AdS boundary. We derive the effective action to check that these 
scalar fields have a well behaved kinetic term. In other words, they are not ghosts. In 
section ini we go to quadratic order in the brane bending. We discover the new scale 
p^, and are able to correct the tensor structure of the propagator. In section [7J we 
briefly show how the mixed case admits self accelerating solutions. Finally, section |S1 




(3) 



3 



contains some discussion. In particular, we argue that the quadratic analysis is valid 
up to very high scales, and that there is no strong coupling problem. 



2 The set-up 

Consider two 5 dimensional spacetimes, A4l and J^r, separated by a domain wall. 
The domain wall is a 3-brane corresponding to our universe. Our set-up is described 
by the following action: 

•S* = 5'bulk + 5'brane (4) 

where the contribution from the bulk is given by 

J2 f d^x^/^ {R - 2Ai) + 2Ml j d^xv^ K^^ (5) 

_r p J Mi J brane 



'bulk 



Here, Mi is the 5-dimensional Planck mass in M.i. We have also included a negative 
cosmological constant, Aj = — 6//f, We will not assume that there is Z2 symmetry 
across the brane, so that the Mj and li can differ on either side of the brane. The 
bulk metric is given by gab, with corresponding Ricci scalar, R. ■jab is the induced 
metric on the brane. The extrinsic curvature of the brane in A^j is given by 

^1? = 7^7.'V(,n,) (6) 

where is the unit normal to the brane in A^j, pointing out of A^j. For the most 
part we will not bother with the index i when referring to bulk quantities, although 
the reader should understand that they are there. 
The brane part of the action is given by 



•'brane 



d'^x {-ay/^ + Cm) (7) 

brane 



where a is the brane tension and Cm describes any additional matter. 
The bulk equations of motion are given by the Einstein equations 

Rab — -^^dab = —^Qab (8) 

The boundary conditions at the brane are governed by the Israel equations [201 • This 
comes from varying the action with respect to the brane metric. Given a quantity 
Zi defined in Jvii we shall henceforth write {Z) = {Zl + Zr)/2, for the average across 
the brane, and AZ = Zl — Zr, for the difference. The brane equations of motion are 
given by 

a 1 /„ 1 



where 



2{M'Kab) - -lab = - {Tab " 'T^ab ) (9) 



2 dC 

Tab = (10) 



is the energy momentum tensor for the additional matter on the brane. 
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3 Background solutions 



In this section we will derive the metric, (jab, for the background spacetime. These 
correspond to solutions to the equations of motion when no additional matter is 
present. Let us introduce coordinates x"" = and assume that the brane is 

located at z = 0. The left hand side, Aii corresponds to z < 0, whereas the right 
hand side, Ain corresponds to z > 0. In order to trust our classical analysis in the 
bulk, we need to assume 

M>1/1 (11) 

Now seek solutions of the form 

ds^ = gabdx^dx^ = a^{z)ri^^dx''dx'' + dz^ (12) 

where rj^^ is four-dimensional Minkowski space. The Einstein equations yield the 
following 

'a'\^ 1 a" _l 

/2' T ~ /2 



(13) 



Without loss of generality we impose the condition a(0) = 1 so that 



a{z) = exp (^y^ (14) 

where = ±1 

Given the solution ()14|) . the Israel equations Q impose the following condition 
on the brane tension 



Note that this gives the usual fine-tuning of the brane tension for the Randall- 
Sundrum model l24i, where AM = A/ = 0, and 6l = ^ = —Or- 



a 



4 Linearised perturbations 

We shall now consider metric perturbations about the background solutions we have 
just derived. We will allow additional matter. Tab, to be present on the brane, but not 
in the bulk. We begin by deriving the bulk equations of motion, and the boundary 
conditions at the brane. 

4.1 Bulk equations of motion and boundary conditions 

Let us define gab = Qab + ^Qab to be the perturbed metric. We will work in Gaussian 
normal (GN) coordinates, so that 

^9tiz = ^Qzz = (16) 
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Since we have no additional bulk matter, we can take the metric to be transverse- 
tracefree in the bulk. In other words, 6g^i, = Xfiu{x,z), where 

d.x; = = x', (17) 

Here indices are raised and lowered with g^^. In this choice of gauge, the linearised 
bulk equations of motion are given by 



52 ^ _ 4 

q2 ^2;2 /2 



X^,uix,z) = 



Unfortunately, we no longer expect the brane to be at 2; = 0. The presence of matter 
causes the brane to bend IHOj- In fact, we would even expect the brane to be bent 
by different amounts when viewed on different sides of the brane. In other words, 
the brane position is given by z = fi{x) when viewed in Aii. This makes it very 
difficult to apply the Israel equations ©. To get round this, we need to apply a 
gauge transformation in Aii and another in Ai^. The new coordinates should also 
be GN, with the brane positioned at z = 0. 

To this end, we make the following coordinate transformations 

nj 

z^z-f{x), x^^x^-g'^(x) + -(l-a-2)5^/ (19) 

Note that in addition to 6 and /, f{x) and Q^{x) are understood to have an invisible 
index i = L,R. In these new coordinates, 

26 

SgtMu = Xi^Ax, z) + 2a^(5(^,^) + 61/ (1 - a^) d^dyf + —fg^^u (20) 

We can evaluate this at z = to derive the brane metric 

16 

h^Lu = XfiA^^ 0) + 2Q(^,^) + —fr]f,y (21) 

For this to be well defined, we require that A (S'jf^y) = 0. Since the geometry is inde- 
pendant of the pure gauge term, 2(5 (^ j,), we should also demand that the remaining 
part of is well defined. Making use of the fact that x'^ = 0, this implies the 
following 

Ax^.u{x, 0) = 0, AQ(^,,) = 0, A (^^^ = (22) 

The last condition suggests that we introduce the function F = 6f /I. The bulk metric 
is now given by 

SgfMu = = Xi^u{,x, z) + 2a^(5(^,^) + (1 - a^) d^,dyF + 2F^^^ (23) 
and the brane metric by 

h,. = hll = xll + 2Qm + 2Fr],y (24) 
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where xjjl/ = Xfiu{x,0). Making use of equation (fT3j). the Israel equations Q imply 
the following boundary condition 

A [M\;,(x, 0)] - ^xll = (25) 

where 

= T^, - ^Tt]^, + 2AiM'ei)d^d,F (26) 
Given that Xfiu is transverse-tracefree, we immediately see that = 0. This gives 

T 

2A{M^ei)d'^F = - (27) 
3 

Here we see explicitly how matter on the brane causes it to bend. 

4.2 Solving the equations of motion 

In four dimensional momentum space, the bulk equation ()18|) becomes 



p2 5^ _ 4 



a 



X>.u{p,z) = (28) 



where tilde represents the Fourier transform. The general solution is 

X,Ap, ^) = A,Ap)h + B,^ip)K2 (29) 

where /„ and Kn are modified Bessel functions of order n |3T] . Naturally, we require 
that Xfj.u{p,z) as \z\ — > oo. To impose this condition we need to consider three 
separate cases: (i) the Randall- Sundrum case for which 6^ = 1 = —Oji, (ii) the inverse 
Randall- Sundrum case for which 6fi = 1 = —6l and (iii) the mixed case for which 
6l = 1 = Or- In figure n we can see the generic behaviour of the warp factor in each 
case. Making use of the fact that [HI] 

l2{y)^0, K2{y)^oo asy^O (30) 
l2{y) oo, K2{y) ^0 as y ^ oo (31) 



we see that the boundary condition at infinity is satisfied if and only if 

-2{pl/a) ~br 
C2{pl) ^'^^ 



X,Ap,z) = ^0^&) (32) 



where 



C2={^' for RS case ^^^^ 
I2 for inverse RS case 



and for the mixed case Cn = Kn and = lo 
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z 

(c) The mixed case. 



Figure 1: The behaviour of the warp factor. For the RS case, a{z) vanishes at infinity. 
For the inverse RS case, a{z) grows exponentially at infinity. For the mixed case, a{z) 
vanishes at — oo but grows exponentially at +00. 
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The Fourier transformed boundary conditions ()25p and ()27|) are given by 



a 



2 ^fiu 



f 
~3" 



Given that [SI] 

yl'^{y) + 2/2(2/) = 2/^2(2/) + 2^^2(2/) 

we can insert the solution (jHSl) into (jSlI) to give 



-yKi{y) 



~br 



R{p) 



where 



Here 



R{p) = -pA 



C2{pl) 



—Ki for RS case 
Ji for inverse RS case 



(34) 
(35) 

(36) 

(37) 
(38) 
(39) 



and for the mixed case Cf = —Ki and = Ji. Note that we have also used equation 
(fT3j) in deriving equation (jHTj) . 

In momentum space, the metric on the brane is given by 



^br ,~br 



(40) 



We can choose Q^{p) to cancel off the term proportional to p^PuF in T,^u{p). This 
leaves 



h 



br 



where 



R{p) ^'^^ 3^''^'^ 



1 



(41) 
(42) 



This all makes sense as long as R{p) and a do not vanish in p > 0. 

At some energy scale, the solution (PT|) should agree with the corresponding result 
from 4-dimensional Einstein gravity 



h 



m'^iP'^ 



T, 



(43) 



This requires there to be an energy scale for which R{p) ~ ap^ and a = m^i > 0. We 
shall now investigate whether or not such a scale exists. 

At energies that are small compared with the bulk curvature, we have ^ 1. In 
this limit we can approximate the Bessel functions as follows |31j 



h{pl) ~ -pi, h{pl) ~ -{pl)\ K,{pl) ~ -, 



K2{pl) 



{plf 



(44) 



At high energies, pi ^ 1, we have 

hjpl) K.jpl) 
h{pl) ^ ' K2{pl) 

For the RS case {6l = 1 = —Or) we see that 



(45) 



fap^ for p < 1//^, 

R{p)=2p(^M^j^ \ ^llMlhp' + Mj^ for 1//^ « p « (46) 

[2(M3)p forp>l//i 



where a = {M^l). Note that we have assumed Ir ^ II- 

At first glance, its appears that we only have Einstein gravity when p <C I/Ir- 
However, if MIIl > MIIr, R{p) ~ \MIIlP^ for 1/Ir -C p < Furthermore, 
we now have a ~ |M|/i. This means that we also have Einstein gravity whenever 
^ P ^ 1//^. In any case, there is no modified gravity in the infra-red. 

For the inverse RS case {Or = 1 = — ^^/^), we see that 



3(MV0 forp«l//« 
R{p)=2p(^M^j^)^\Mlp + A^ forl//^«p«l//i (47) 
^ 2{M^)p forp>l//L 



and a = —{M^l). Einstein gravity cannot be produced at any scale. 
Finally, for the mixed case {Ol = 1 = Or), we see that 



R{p) = p 



IMIIlp" + 4^ for p < 1/Ir 
\MIIlp'' + MIp for 1/Ir < p < l/h 
2{M^)p forp>l//L 

(48) 

and a = ^A{M^l). Again, at first glance its appears that Einstein gravity can never 
be reproduced. However, as in the RS case, if M|/l ^ M^Ir, we have Einstein 
gravity for 1/Ir ^ p ^ The difference here is that gravity can be modified in 

the far infra-red region, p ^ Unfortunately, as we shall see in sectional the 

mixed case contains a radion ghost whenever a > 0. Here a ~ }^M^Il > 0. 

These results are a little disappointing. We can only get Einstein gravity in the 
RS case, without any modifications in the infra-red. This is nothing new [21]. 

However, all is not lost! Consider the RS case and the mixed case for p ^ 1//^- 
Now suppose Mj^li -C M^Ir, in contrast to what we discussed earlier. This gives 

f \MIIr for the RS case , , 

« ~ < 1 , (49) 
I — ^M^/ij for the mixed case 

Since a < for the mixed case, there is no radion ghost. In both cases, we see the 
introduction of a new length scale 



M'r 



< Ir (50) 
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We now have 

'M|j9 for < p < 1/r 



f 

R{p) ~ <^ \MIIlp^ for 1/r < p < 1//^ (51) 
[m|j9 forp>l//L 



Note that we have also assumed that ^ M^. In these hmits, the brane bending 

''fii/- 



term, T /ap^ is always much less than the remaining terms in h!°J^. This implies that 



/i*^' ^ T - -Tn 

11^,,,, ~ „/ N ^1/ „ -t '/^jy 



(52) 



In the intermediate regime (l/r<^p<^l//i), our solution has the correct momentum 
dependence, but the wrong tensor structure. As we move into the infra-red ^ 
p ^ 1/r), or the ultra-violet {p ^ 1/^l) our momentum dependence changes so that 
gravity is modified. 

Of course, our solution (j52|) no longer makes sense when /ij^^ is of order one. For 
a source of mass, m, note that |T| ~ mp^. This means that /ijj^ schematically goes 
like mp^/R{p). This becomes of order one when p ~ Pcut-ofr- For very large masses 
(m > MflLf), this cut-off occurs in the infra-red. Otherwise we have 

_ {Mlh/m for Mill < m< Mlhr 

Note that the upper value in (jSHj) lies in the intermediate energy range, whereas the 
lower value lies in the ultra-violet. 

We might think that this cut-off represents the scale at which the linearised theory 
breaks down. Certainly, the linearised theory makes no sense beyond the cut-off, but 
is it really valid beforehand? In the DGP model the linearised theory actually 
breaks down much sooner than expected. This is because higher order terms become 
important ^H] • Curiously, these non-linear terms can sometimes correct the problems 
with the tensor structure at intermediate energies ISHl IMl tHHl ISE] ■ In section El 
we will see that the same thing happens here. 



5 Scalar and vector perturbations 

We shall now discuss scalar and vector perturbations in the bulk, when there is no 
additional matter on the brane {Tab = 0). 

5.1 Scalar perturbations and ghosts 

Braneworld models that exhibit infra-red modifications of gravity often contain ghosts 
(see, for example |lll2ni)- Typically, these ghosts appear in the scalar sector. We will 
be forced to eliminate any model that contains a ghost-like scalar field in the effective 
theory. 
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Consider the following scalar perturbations about the background, gab, 



6g^, = a" {2d^d,E + 2Ar,^,) , 
Under the scalar gauge transformations 

x'^ + 

the scalar fields transform as follows 

e 



^9, 



6gzz = 20 



(54) 



(55) 



i?-(e. + aT), (56) 



It is easy to check that the following quantities are gauge invariant 

X = - ieA\ Y = B- a^E' - WA 



(57) 



The boundary condition at infinity requires there to be no physical perturbation there. 
This amounts to the gauge invariants X and Y vanishing at infinity. We can see this 
most easily by choosing the gauge E = A = so that X = cf) and Y = B. 

In terms of the gauge invariants, the bulk equations of motion (jHl) are given by 



= 
= 
= 

The solution is 



36 



d^Y' 



Y' 
d„X 



29, 



-Y -2X 



9 



a'' I 



I 



X = 0, 



/2 



Y 



U{x) 



(58) 
(59) 
(60) 

(61) 



where d^U = 0. In principle, U can differ on either side of the brane. 

Recall that we require X, F — > as \z\ oo. For the RS case, this requires that 
Ul{x) = Ur{x) = 0. For the inverse RS case, U{x) can be non-zero on both sides of 
the brane. Finally, for the mixed case, we must have Ul{x) = 0, although Uji{x) can 
be non-zero. 

In order to apply the boundary conditions near the brane we have to choose a 
gauge. We can choose GN gauge whilst keeping the brane position fixed at 2; = 0. 
We now have B = (j) = 0, and 



i^/f/(x)(a-^ - 1) + h^V{x){a'^ 



E = -9lU(x)(a-'' -1) + =:FV(x)(a'' -1) + W{x), A = V{x) (62) 
We can evaluate our solution at 2; = to derive the brane metric 

57/.. = 2d^d,W + 2Vri^, (63) 

The pure gauge part, 2d^duW ^ and the remainder, 2Vrj^y, should both be well defined. 
This means that AV^ = /^W = 0. In particular, since AW = 0, it is easy to see that 
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W can be continuously gauged away on both sides of the brane. We therefore set 
W = 0. 

The Israel equations Q are now given by 

2df,d^A{M^E') = (64) 

where we have made use of equation pSj) . This implies that 

A{M^U) = -A{M^ei)V = -2aV (65) 

For the RS case {9l = 1 = —Or), we have deduced that Ul = Ur = 0. Since 
a = (M^l) 7^ 0, equation ()65|) clearly implies that = 0. This means that there is 
no scalar perturbation in the RS case. 

In contrast, for the mixed case and the inverse RS case, we conclude that there 
is at least one remaining degree of freedom, V, say. This means that there is at least 
one scalar perturbation, which we will refer to as the radion. 

We need to check whether or not the radion is a ghost. This involves calculating 
the radion effective action. Again, we need to choose a gauge. Our gauge choice 
must correspond to the brane being fixed at z = 0, and the perturbation vanishing 
at infinity. 

Let us begin by choosing the gauge E = A = 0, for which 

B = ^, = (66) 

This vanishes at infinity. Near the brane, we begin in GN gauge, but transform to 
E = A = hj choosing 

e = ^eiU{x)ia-^ - 1) + ^/V(x)(a-2 - 1), = lOV (67) 

The brane is now positioned at z = WV . We need to move it back to z = Q without 
introducing a perturbation at infinity. This can be done by choosing ^2 = P{z), where 
P(0) = —19V and P ^ {) as \z\ ^ 00. The bulk perturbation is now given by the 
following 

I a 
To quadratic order, the effective action is 



E = 0, A = --P{z), B = --P{z), <P=-P\z) (6^ 



-^eff =-\ j V^M'h''''5Eat - \ f V^h'^'Se^t, (69) 

^ ^bulk ^ J brane 

where 6Eab and 6Qab are the expansions, to linear order, of the bulk and boundary 
equations of motion respectively 

Eab = Rab - ^Rgab + J^Qab (70) 
Qab = 2{M'{Kab-K^ab)) + ^7ab (71) 
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We find that 
5'efr = 



-? [ d'x ^d'U, r dzP'iz) + ^d'Un rdzP\z) 

J II J-oo Jo+ 

^ j d^xVd^A{M^U) = -3a j d^xVd^V 



(72) 



where we have made use of equation We immediately see that the radion is 

a ghost whenever a > 0. This is not a problem for the inverse RS case, as a = 
— {M^l) < 0. However, as we saw in the last section, this can be a problem for the 
mixed case. 



5.2 Vector perturbations 

Now consider vector perturbations in the bulk. Since brane bending represents a 
scalar fluctuation we can set the brane to be at 2 = 0. In general, the vector pertur- 
bation is given by 

6gf,„ = 2a^F(^,^), = 6g^^ = (73) 

where and -B^ are divergence-free. 

Under the vector gauge transformation, 

^xf' + C, z-^z (74) 
the vector fields transform as follows 

F^^F^.- C,, B,^B,- aX (75) 

where = t]^^^^ . Note that is divergence-free. It is clear that we have the 
following gauge invariant 

C, = B,- a'Fl (76) 
The bulk equations of motion (jH} are given by 

= + (77) 



2^ 



= -o-C/. (7^ 



These have the following solution 



C, = ^ (79) 



where c^^A^ = 0. Of course can differ on either side of the brane. In analogy with 
the scalar gauge invariant, we require that — > as \z\ — > 00. For the RS case, this 
requires that A^ = A^ = 0. For the inverse RS case, A^ can be non-zero on both sides 
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of the brane. Finally, for the mixed case, we must have = 0, although can, in 
principle, be non-zero. 

To apply the boundary conditions we go to GN gauge (5^ = 0), so that 

F^ = -^ei\^{a-^~l) + G^{x) (80) 

Continuity at the brane implies that AG^ = 0, so we can continuously gauge away 
Gfj,. The Israel equations give 

2A [M^F^,,)] = (81) 

This implies that 

A [M^A^] = (82) 

It is easy to see that if A^ is zero anywhere, it is zero ef er?/where. We conclude that 
vector perturbations can only exist for the inverse RS case. 



6 Beyond the linearised analysis 

At the end of section HI we entertained the possibility of the linearised analysis break- 
ing down at an unexpectedly low scale. In this section, we will show that this is 
indeed the case. We will go to quadratic order to gain insight into the full non-linear 
theory. 

Let us begin with the transverse-tracefree + GN gauge in the bulk 

= Xiiu, Sg^^ = 5g^^ = (83) 

The brane is now bent. We need to fix the position of the brane to be at 2 = 0, 
whilst remaining in GN gauge to quadratic order. To this end, we make the following 
coordinate transformation 

z^z-f{x)--{l-a-^)d,fd^f (84) 



9/9 

d^f + dJd^Q'' + -fd^f 



+ 



X^^'dJ (85) 



Note that this agrees with the linearised transformations (fT!I|l to first order. The bulk 
metric is now given by 

5g^u = Ku = Xt^Ax, z)+2a2Q(^,,)+/2 (1 _ d^d,F+2Fg^,-Pdf,Fd,F+6h^, (86) 

where F = Of /I. The terms in 5h^y are always much smaller than the other terms 
in h^ij^ so we will neglect them. For the linearised analysis to be valid, we require 
that the linear terms in ()86|1 are much larger than the remaining quadratic term, 
l^d^Fd.F. 
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Suppose p ^ 1/^_R, Ml S> Mr and MfZ^ <^ Mf^l^. For a source of mass m < 
Mfl^r, we have shown that the hnearised analysis certainly breaks down when p ~ 
Pcut-off, where Pcut-ofr is given by equation We will now show that the quadratic 
term, t^d^Fd^F, actually becomes important much sooner. 

From sectional we know that the linear terms go like mp^/R{p). Schematically, 
the quadratic term becomes important when p^t^F"^ ~ mp^ /R{p). To identify when 
this happens we need to know the size of F. From equation (j35p . we see that 

1^1 1^1 (87) 



' ' miirp^ miir 

It turns out that the quadratic term becomes important when p p^ <ti Pcut-os- 
For the remainder of this section we will restrict our attention to sources with the 
following mass 

Mlll/r^ <m< MIIlt (88) 

When r is large, this covers a huge range of masses. By taking our mass to lie in this 
range we ensure that 1/r < p^ < Specifically, 

v.- mi m 

We conclude that the linearised analysis only makes sense for p <^ p^ ^ Pcut-ofr- For 
p* <^ p <^ Pcut-off, we cannot ignore the quadratic term, Ij^d^F^duFR. Let us proceed 
with this in mind. 

It is convenient to decompose the remaining quadratic as follows 

-fd^Fd.F = T^, + 2F(^.,,) + Id^d.E + 2Ar/^, (90) 

where JF^jy is transverse-tracefree and is divergence free. It is easy to show that 

d^A = — Id^d" - r]^''d^] d.Fd.F (91) 
6 

Since jF^i, is transverse-tracefree, let us absorb it into a redefinition of X/iu- Further- 
more, let 

^ - - d^.E (92) 

so that the bulk metric is given by 

V = Xt^uix, z) + 2a2Q(^^,) + /2 (1 _ a?-) d^d^F + 2Fg^, 

+ 2(1 - a^) + d^d^E] + 2Ar]^, (93) 

We can evaluate this at 2; = to derive the brane metric, 

^7m. = = x^.u{x, 0) + 2g(^,,) + 2(F + A)ri^, (94) 

The pure gauge part, 2(5 (^ j,), and the remainder must be well defined. Making use 
of the fact that = 0, we find 

Ax^.(x,0)=0, AQ(^,,) = 0, A{F + A) = Q (95) 
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It turns out that the Israel equations take very nearly the same form as in the lin- 
earised theory, the only subtely being that AF ^ 0. Once again we have equation 
(1^ but with 

= T^, - ^Tr]^, + 2d^d,A{M^eiF) (96) 
Note that we have used the fact that 

Fi^,,) + d^d.E + A'q^, ~ pH^F^ < pH^F ~ fd^d.F (97) 
Given that x^v is transverse-tracefree, we have that 

2d'^A{M^eiF) = - (98) 
3 

In momentum space, the metric on the brane is given by 

hll = + 2*P(M^-) + + A)v,u. (99) 
As in section HI we see that 

xt = %¥ (100) 



We choose Q^{p) to cancel off the term proportional to p^pyA{M^OlF) in S^j,(p). 
This leaves 

1^ 



^"•^ R{p) f 3 
For the RS case and the mixed case, recall that 



+ 2{F + A)r]^, (101) 



m-l^H . (102) 
\\MIIlp^ forl/r<p<p,, 

provided Ml > Mr and M|/i < M^lr. Here 

Puv = min{l//L,pcut-ofr} = min{l//L, MlZi/m}. (103) 

Our goal is to reproduce Einstein gravity at certain scales. R{p) already has the 
correct momentum dependance when 1/r <^ p <^ p^^. To ensure that our solution 
also has the right tensor structure (see equation pH|l ). we require that 

f 

F + A = ^ (104) 

We will now see that this is indeed the case when p^ <^ p <^ p^v 

For p <^ p*, we know that the linearised analysis can be trusted. However, for 

p* <^p <^ Puv, we have 

Fl > Al, Fr < Ar (105) 
Making use of boundary condition (jHH), we deduce that 

Fl^ Fl + Al = Fr + Ar^ Ar (106) 
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Now consider the trace equation ()98|). In principle, we have two possibihties: (i) 
MIIl\Fl\ > MIIr\Fr\ or (ii) MIIl\Fl\ < MIIr\Fr\. For < p < p^,, it turns 
out that case (i) is the only one to give a self-consistent solution. The trace equation 
now implies that 

-m|i? <^°^) 

SO that equation ()l()4j) holds. 

To briefly sum up, we have shown that the RS case and the mixed case mimic 
four- dimensional General Relativity when -C p <^ This happens because the 
linearised results break down, and we have to include quadratic contributions. At 
lower energies, the linearised results hold, and we deviate from four dimensional GR. 
At first, the tensor structure of our propagator changes. As we go even deeper into 
the infra-red, our momentum dependance also changes. 

It is important to ask whether or not the radion analysis of section^lis significantly 
affected by non-linear effects. The answer is clearly "no" . This is because the radion 
field is massless = 0, so the linearised analysis can be trusted. Recall that there 
is no radion in the RS case, and the radion is not a ghost for the mixed case when 
Mf^li < MrIr. We conclude that the radion is not a problem here. 



7 A word on cosmology 

Infra-red modifications of gravity generically have far reaching implications for cos- 
mology, particularly at late times. In this section, we will briefly examine the cosmo- 
logical solutions that arise in our models. We will omit most of the details, as they 
can be found in We should note that in [21], Or is defined with the opposite sign 
to here. 

Let us consider a spatially flat Friedmann- Robert son- Walker brane sandwiched in 
between two copies of anti-de Sitter space. The dynamics of the brane are governed 
by the following equation [2lj 



A 



^ (108) 



where H is the Hubble parameter on the brane, and p is the energy density. Here p 
is made up of brane tension and additional matter 

p = cr + p^ (109) 

Note that a is given by equation (fT3j). For 1/Ir_ ^ H <^ 1/^l, equation ()108j) 
approximates to 

OlH^ - 29 J-^ ^ (110) 



At early times, H ^ 1/r, so that 

^2^0LPrn ^ ^ 

3MflL ^ ^ 
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For the RS and the mixed case, 6*^ = 1, so this simply gives the standard 4d cosmology 
with mil = \MIIl. 

We might expect to see a de Sitter phase of expansion at later times, as — ^ 0. 
For the RS case {6l = —Or = 1), it is clear from equation (jlOSj) that this will not 
happen: = is the only solution for pm = 0. In contrast, for the mixed case 
{9l = 9ji = 1), the solution H ^ 2/r is permitted. 

We conclude that late time de Sitter expansion occurs for the mixed case, but not 
the RS case. The RS result is a surprise as we have modified gravity in the range 
1/Ir < p < 1/r. 

8 Discussion 

We have considered in detail perturbations about a fiat brane embedded in between 
two copies of anti-de Sitter space. We have abandoned Z2 symmetry across the brane, 
giving rise to three distinct cases: (i) the RS case where the warp factor decays away 
from the brane on both sides, (ii) the inverse RS case where the warp factor grows 
on both sides and (iii) the mixed case where the warp factor decays to the left of the 
brane, and grows to the right. 

Let us now summarize our main results. The inverse RS case is a dead end: four 
dimensional Einstein gravity can never be achieved. In contrast, take the RS case 
and the mixed case when S> Mji and M^l^ <C Mf^l^. Consider a source of mass, 
m G {Mfl\/r'^^ M^l^r), where r = M^li/M^. The metric on the brane is given by 

when < p < Puv 
when 1/r <^ p <^ (112) 
when <^ p <^ \/r 

where m^, = ^MIIl, p* = {MjlL/mr'^Y^'^ , and puv = minjl/Zi, M|/i/m}. Gravity 
on the brane mimics four-dimensional GR at intermediate scales, but is modified 
in the infra-red. This modification begins with a change in the tensor structure of 
the propagator. At even lower energies the momentum dependance also changes. 
Crucially, these models are free from ghosts. 

We should also note that a cosmological brane in the mixed case will approach a 
de Sitter phase at late times, even when there is no additional matter. This does not 
happen for the RS case. 

To illustrate our results more clearly, let us put some numbers in. Consider 
the mixed case with Ml ~ l/h ~ rupi ~ 10^^ GeV, Mr ~ 10 - 100 MeV, and 
1/Ir ^ 10"^"^ eV. Firstly, note that the current Hubble rate, Hq ^ 2/r ~ lO"^'' eV, 
as desired. Now consider our source to be the Sun, with mass m ~ 10^^ eV. Four- 
dimensional GR is reproduced at distances between rschw and 10^^ cm. Here, the 
Schwarzschild radius of the Sun is given by rschw ~ 1 km. The large distance limit, 
10^^ cm, extends well beyond the size of the solar system, so we have no conflict with 
experiment. 



''fii/ 



T, 
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Of course, the results ()112|) were derived by going to quadratic order in perturba- 
tion theory. When <^ p <^ Puv, we found that the quadratic brane bending term 
dominated over the hnear terms. We might be worried that this means we have lost 
control of perturbation theory. However, we will now argue that we do have control 
of perturbation theory, because terms beyond quadratic order are always small, right 
up to Puv 

Suppose, we have two GN coordinate systems, x"" and where 

= x^' + i^{x,z), z = z + C{x,z) (113) 

In the "unhatted" coordinates, the brane is positioned at z = 0. For simplicity we 
will demand that C,^{x, 0) = 0. This just corresponds to a choice of gauge on the 
brane. In the "hatted" coordinates the brane is positioned at z = ^^(x, 0) = f{x). 
Non-perturbatively, the metrics in each coordinate system are related as follows 

Ox^ dx'^ 

Qabix, z) = gcd{x, ^)-Q^-Q^ (114) 
We shall now evaluate the (/ii/) component at 2; = to derive the metric on the brane 

Qi^uix, 0) = ^^^.(x, /) + d^fdj (115) 

Suppose we write g^y = g^^ + x^lV1 and Taylor expand the right-hand side of (|115|) . 

We get g^,y{x, 0) = g'l^'"^{x) + 6g^^{x), where 

gTu'^ix) = V(x,0) + /^;,(x,0) + ^/2^;:,(a;,0)+x^.(x,0) + /x;,,(a:,0)+9^/9,/ (116) 
and 

°° -Cn ^ -Cn 

6g,.{x) = J2 ^ai-Jix, 0) + J2 0) (117) 

n=3 n=2 

Here prime denotes partial differentiation with respect to z. g'j[^^'^{x) includes all those 
terms that were considered in the analysis to quadratic order. Sgfj,u{x) contains all 
those terms that were assumed to be of higher order, and were consequently ignored. 
We shall now ask whether or not it was safe to ignore them. 

Roughly speaking, we can associate each z derivative with either a factor of p or a 
factor of 1//. It is obvious that as long as p/ <C 1 and /// -C 1, then Sg^u <^ g'ji'if''^- In 
the quadratic analysis it is clear that pf and f /I are indeed small whenever p -C Puv 
We conclude that it was entirely consistent to ignore all terms beyond quadratic order. 

In analogy with the DGP model, we shall assume that a break down in pertur- 
bation theory gives rise to strong coupling. To derive the strong coupling scale we 
consider a source of mass m ~ rUpi, and ask at what scale do we lose perturbative 
control j2Hl- If this scale were p^, we would be in trouble. We would expect strong 
coupling below distances of around 1000 km. However, we have just shown that we 
don't expect to lose control until at least Puv, which typically corresponds to the 
Planck scale! These results suggest that there is no strong coupling problem here, 
although a more thorough analysis is clearly required. 
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Finally, note that the mixed case in particular appears to have much in common 
with the DGP model The expressions for r and p^, seem to have obvious DGP 
analogues. The cosmological behaviour is also very similar (see equation (jliup ). One 
can't help thinking that this is more than mere coincidence. With this in mind, con- 
sider the brane to be the common boundary to the decaying bulk and to the growing 
bulk. If we were to calculate the boundary stress-energy tensor on the decaying side, 
we would expect there to be divergences. These must be cancelled by boundary 
counter terms that will probably take the form of localised curvature on the brane. 
If this is indeed the case, it would represent a new mechanism for obtaining induced 
curvature on a DGP brane. 
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